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1. Introduction 

The method of large N continuum reduction [J], for SU (N) gauge theory allows for the 
calculation of the infinite volume, infinite N limit of certain physical quantities using volumes 
reduced to a small physical size. Numerical estimates [JjJ ||] of the physical critical size above 
which continuum reduction holds indicate that this method can be used to produce practical results. 
The chiral condensate ^ and the pion decay constant[|j] were calculated in the large N limit 
in four dimensions using continuum reduction. In this paper, we show that the method can be 
extended beyond bulk quantities and that it also produces reliable results for quantities with space- 
time dependence such as the heavy quark potential, from which the string tension can be extracted. 
Invoking large N continuum reduction, we included Wilson loops larger than the size of the lattice. 
The results validate the method of continuum reduction for calculating quantities based on the 
space-time dependence Wilson loops. 

A precise calculation of the string tension in three dimensional SU (N) gauge theories has 
been performed with N up to 8 on large lattices [Q]. We present a complementary calculation with 
N = 47 on 5 3 lattices using continuum reduction. The calculation of Ref. ^ used correlation 
functions of smeared Polyakov loops to extract the string tension. After extrapolating to N = °° and 
to the continuum, the result was 

0.1975 ± 0.0002 - 0.0005 (1.1) 



8 2 N 



where g is the gauge coupling. This has to be compared with the analytical calculation in [|7|], 
namely, -4= « 0.1995. Although the two results are not in perfect agreement, the main observation 



is that the approximations used in the analytical calculation are very well motivated. 

Our use of continuum reduction to directly compute the N = 00 limit of the string tension by 
working at large enough N so that the finite N corrections are smaller than the numerical errors 
gives [@| 

Xr- =0.1964 ± 0.0009. (1.2) 

g N 

This result and that of ( |IT| ) are consistent at the level of their one sigma errors. This level of 
agreement is, in turn, consistent with neither the large N extrapolation of Ref. []|] nor the volume 
reduction of the present calculation having unexpected errors. While both of the numerical results 
lie below the analytical estimate, the discrepancy is relatively small. Thus the numerical evidence 
that the analytical result is an excellent first approximation that captures much of the physics re- 
mains strong. 

The paper is organized as follows. We explain how we use smeared Wilson loops to compute 
the string tension in Section ||. The lattice results for the string tension along with the continuum 
extrapolation are also presented in this section. An intermediate step in our calculation is the 
dimensionless ground state string energy m{k). In Section ||, we show results for m(k) at one 
fixed lattice coupling to illustrate its behavior as a function of k and how it is used to extract the 
string tension. We also show that m(k) is unaffected by the smearing parameter. We illustrate the 
extraction of m(k) at one fixed coupling in Section ^ Here we show how the smearing parameter 
affects the overlap with the ground state. The main result in this paper is obtained using N = 47. 
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We show that the finite N and finite volume corrections are small at this value of N in Section ||. 
We explain why this method is preferred over the Creutz ratio in Section |6| 

2. String tension using Wilson loops and continuum reduction 

Consider SU (N) Yang-Mills theory on a periodic lattice with the standard Wilson gauge ac- 
tion. The method of [Q] is to measure the string tension using correlations of Polyakov loops with 
separation t that wind around a space direction. Continuum reduction [jlj |2|] implies that the large 
N Yang-Mills theory in a continuum box of size Z 3 is independent of / as long as / > l c = 1 /T c with 
T c being the deconfining temperature. One should be able to compute expectation values of Wilson 
loops of arbitrary size on an I 3 continuum box using folded Wilson loops and extract the string 
tension. To implement this approach to the three-dimensional Yang-Mills theory string tension, we 
use the following procedure: 

• We fix the lattice size to L 3 . We use L = 5 for the most part and only use L = 4 to verify 
reduction. 

• We fix N so that finite N corrections are small. We set N = 41 and show using one instance 
that finite N corrections are small at N = 47. 

• We pick an appropriate range of lattice coupling b = 

- b cannot be too small since we have to be away from the bulk transition on the lattice 
associated with the development of gap in the eigenvalue distribution of the plaquette 
operator [§]. Therefore, we pick b > 0.6. 

- b cannot be too big since we have to be below the deconfining transition for L = 5. 
Therefore, we pick b < 0.8 [§]. 

• We use smeared space-like links and unsmeared time-like links. 

• We use the tadpole improved coupling bj = be(b) to set the scale and consider K xT Wilson 
loops W(K,T) with 1.5 < ^, Jj < 12.5. This amounts to expectation values of Wilson loops 
that range from 0.82 to 2 • 10~ 4 . 

• Keeping K fixed, we fit 

lnW(k,t) = —a — m(k)t; (2.1) 

where k = and t = Jj are the dimensionless extent in the space and time direction respec- 
tively. m(k) is the dimensionless ground state energy. This fit assumes that there is a perfect 
overlap with the ground state. Note that a should be zero since W(k,0) = 1. Any small 
deviation from zero seen in the fit is due to the contribution from excited states. 

• Finally, m(k) is fit to abjk + cobj + The combination yfcbi is plotted as a function of bj 2 . 
We expect lattice spacing effects to lead off as by 2 in Yang-Mills theories and this is indeed 
the case in Fig. [j] The continuum limit extracted from this figure was quoted in Eqn.( 1 .2). 
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Figure 1: The string tension is plotted as a 
function of the lattice spacing bj 1 . The fit 
is an extrapolation to the continuum. 




k 



Figure 2: The ground state energy m(k) as 
a function of k for the coarse and fine lattice 
spacings considered here. 



The use of smeared links improves the measurement of Wilson loops. They enhance the 
overlap of the space-like sides of the Wilson loops with the ground state. This increases the signal 
relative to the fluctuations and simplifies the t behavior of the loops [ 1C]. One step in the iteration 
takes one from a set Zj[ (x\,X2,t) to a set uf; l+1 \xi,X2,t). Before reunitarization, the weight of 
Zjff* (x\,X2,t) is (1 — /) while that of each staple is f/2. The time-like links, Uj(xi, x 2 ,t), are not 
smeared, and the smearing only involves space-like staples. There are two parameters, namely, the 
smearing factor / and the number of smearing steps n. Only the product z = fn matters, and / 
plays the role of a discrete smearing step. For a given z, the overlap of the smeared loop with the 
ground state does not depend on / as long as it is small. But the overlap of the smeared loop with 
the ground state does depend upon z. We set the value of the smearing parameter to z = 2.5 by 
choosing / = 0. 1 and n = 25. To study the effect of varying z, we also consider z = 1 .25 (f = 0.05 
and n = 25) at one coupling. 



3. Extraction of string tension 

SU (N) gauge fields were generated on a 5 3 periodic lattice using the standard Wilson action. 
One gauge field update of the whole lattice [Q] is one Cabibbo-Marinari heat-bath update of the 
whole lattice followed by one SU (N) over-relaxation update of the whole lattice. A total of 1500 
such updates were used to achieve thermalization. Measurements were separated by 10 such up- 
dates and all estimates are from a total of 832 such measurements. Errors in all quantities at a fixed 
b and N were obtained by jackknife with single elimination. 

The ground state energy m(k) obtained as a function of k = Jp is fit to 

m(k) = abjk + c b } + ^- (3.1) 
k 

We expect abj to approach a finite value in the continuum limit (bj — * °°). The three parameter fit 
of m(k) as a function of k is shown in Fig. |2[ 



4. Extraction of m(k) 

The dimensionless ground state energy m(k) is extracted at a fixed k by fitting \nW(k,t) to 
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Figure 3: Plot of lnW(k,t ) as a function of 
t for seven different values of k at b = 0.8 
with x = 2.5. 



Figure 4: Plot of \nW(k,t) as a function of 
t for seven different values of k at b = 0.8 
with t = 1.25. 



—a — m(k)t as discussed in Sec. ||| While m{k) should be independent of the smearing parameter 
T = fn, the value of a is expected to depend on x. 

We will use b = 0.8 as the coupling to illustrate the extraction of m(k). Figure [| and Fig. ^ 
show the performance of the fit for two different values of X, namely, 2.5 and 1.25 respectively. 
The solid circles show the data points without errors. The solid lines show the fit of the data. Seven 
values of t were used to fit the data at one k, and data at seven different values of k were fitted. This 
amounted to all Wilson loops from 1 x 1 to 7 x 7 on the 5 3 lattice. The fit parameters are shown in 
Table and Table §. Only the average values of the fit parameters are listed. 

Investigation of Table [T] and Table || shows that m(k) does not depend on x. There is a small 
difference in the two values of m(k) at a fixed k for the two different values of x if k is large. 
Additional analysis shows that this difference is within errors. Furthermore, the fitted values of 
abj for the two different values of X are the same within errors. 



k 
a 

m(k) 



1.62 
0.001 
0.133 



3.23 
0.003 
0.218 



4.85 

0.009 

0.286 



6.47 
0.019 
0.347 



8.08 
0.055 
0.399 



9.70 
0.047 
0.464 



11.31 
0.071 
0.517 



Table 1: Fit parameters corresponding to the fit lnW(k,t) = —a — m(k)t for seven different values of k at 
b = 0.8 with T = 2.5. 



k 1.62 3.23 4.85 6.47 8.08 9.70 11.31 
a 0.002 0.012 0.029 0.054 0.102 0.114 0.144 
m(k) 0.133 0.218 0.287 0.349 0.404 0.468 0.526 

Table 2: Fit parameters corresponding to the fit lnW(k,t) = — a — m(k)t for seven different values of k at 
b = 0.8 with T= 1.25. 

The values of a in Table [T] and Table || do show a variation with x and k. Since a smaller 
value of x implies less smearing, the overlap with the ground state is less for smaller X, and this 
results in a larger value of a at smaller x. The value of a is very close to zero for small k indicating 
excellent overlap with the ground state for the chosen value of x. As k increases, the length of the 
loop increases and the perimeter divergence has a stronger effect. This results in a larger value of 
a as k increases at a fixed x. 
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5. Finite N effects 



Two issues need to be addressed with the analysis performed so far. We have fixed our value 
of N assuming finite N effects are small. If N is not large enough, finite N effects need to be 
addressed. In addition, we also have to address finite volume effects since continuum reduction is 
valid only in the N — > °° limit. 

We expect m(k) to have a fixed limit as N — > °° at a fixed k, L, b and z. Indeed, this is the 
case as shown in Fig. ^ where the results for m{k) as a function of k are shown for b = 0.8 with 
T = 2.5 on 5 3 lattice. All three fit parameters are consistent within errors all the way from N = 23 
to N = 47. The only glitch one sees is at k 8. This corresponds to K = kbj = 5, which is the 
linear extent of the lattice. One can argue that there are larger finite N effects at strong coupling for 
K = L. Since the fit of m(k) involves several values of k, the larger effect at this particular value of 
k is diminished in the extraction of abj. 

Since finite N effects can be ignored at N = 47, we also expect there to be no appreciable finite 
volume effects at this value of N. This point is illustrated in Fig. |6] where the result for m(k) is 
plotted at Z? = 0.6 and % = 2.5 on 4 3 and 5 3 lattice. We used b = 0.6 for this comparison since we 
have to be in the confined phase on 4 3 lattice. Figure ^ shows that the two values of m(k) at a fixed 
k are consistent with each other within errors. The same is the case for the fit parameter abj. This 
is not the case for c\ and c$bj, and this is probably due to a three parameter fit using only five data 
points. Sub-leading coefficients are expected to depend sensitively on the data points. Since we are 
primarily concerned with the value of the string tension in this paper and since all our results are 
based on data taken on 5 3 , we expect the final result to be free of finite Af and finite L errors. 



6. Creutz ratio 



It is natural to ask how the Creutz ratio wM, 



/tV ' ' W{K,J-\)W{K-\,jy 

performs as an observable from which to extract the string tension. If we were to use Creutz ratios, 
we would have smeared all links using all staples. But one can still ask how the Creutz ratio behaves 
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with the asymmetrically smeared links. The K x K square Creutz ratios do not converge well as K 
increases. It is possible the situation would be different if we had smeared all links. 

Each data point in a Creutz ratio is obtained using only four different Wilson loops, i.e. four 
of the data points in Fig. ||. This is quite different from the analysis in this paper. Seven different 
Wilson loops in Fig. || are used to extract one m[k) point in Fig. |2[ and the loops used for different 
k form independent sets. Then the m[k) are fit to determine the string tension. Both folded and 
unfolded loops contribute together. This is the main reason we succeeded in extracting the string 
tension using the range of Wilson loops considered here. To extract the string tension using Creutz 
ratios, larger loops and therefore larger statistics and possibly larger N would be needed. 
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